A Model for Bone Strength and Osteoporotic Fractures 
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Inner porous regions play a critical role in the load bearing capability of large bones. We show that 
an extension of disordered elastic networks [Chung et. al., Phys. Rev. B, 54, 15094 (1996)] exhibits 
analogs of several known mechanical features of bone. The "stress-backbones" and histograms of 
stress distributions for healthy and weak networks are shown to be qualitatively different. A hereto 
untested relationship between bone density and bone strength is presented. 



Osteoporosis is a multi-faceted metabolic disease that 
reduces bone strength and leads to a significant num- 
ber of fractures occurring in older adults |l|,|[. Unfor- 
tunately, therapeutic agents available for prevention and 
treatment of osteoporosis often induce adverse effects in 
patients Q . Thus, non-invasive diagnostic tools to deter- 
mine the need for therapeutic intervention are essential 
for effective management of osteoporosis. Simple models 
can form a useful complement to traditional studies of 
osteoporosis. In this Letter we introduce such a model, 
and describe some of its features. 

Large bones such as thigh bones and vertebrae consist 
of an outer cylindrical shaft (cortex) and an inner porous 
region (trabecular architecture) g]. The cortex is made 
of compact bone and has a thickness of several millime- 
ters. The structure of the trabecular architecture (TA) is 
that of a disordered cubic network of "trabeculae" whose 
axial and cross sectional dimensions are of the order 1mm 
and 0.1mm respectively H (see Figure 1). 

Routine activities (e.g., climbing stairs) inflict micro- 
damage on bone. Material in the neighborhood of these 
fractures is resorbed by a class of cells known as "osteo- 
clasts." Their presence attracts a second group of cells, 
"osteoblasts," which help regenerate lost bone [0. This 
sequence of events, referred to as bone remodeling, re- 
duces the accumulation of micro-damage shallower than 
about 0.1mm p[0|. The full restoration of bone strength 
typically takes a period of 2-3 months. 

In the cortex, deeper fractures, created during occa- 
sional trauma, are not repaired through remodeling. The 
resulting micro- fracture accumulation leads to lower bone 
quality and fracture toughness P,p|, reducing the load 
bearing capacity of the cortex with aging. 

Since the thickness of trabecalae is ~0.1 mm, fractures 
that do not sever them can be remodeled. Clinical stud- 
ies also indicate that perforated trabeculae are seldom 
regenerated B and that cross sections of those surviving 
changes little M. Thus, even though the connectivity of 
the TA reduces with aging, the surviving trabeculae can 



(mostly) be expected to retain their quality. Indirect evi- 
dence for this conjecture is provided by mechanical stud- 
ies which have shown that a TA from a given skeletal 
location fractures at a fixed level of strain independent of 
the age of the bone S. (In contrast, the breaking stresses 
of such samples reduce with age, though much less than 
the corresponding degradation of a cortex.) Due to ef- 
ficacy of its remodeling, the TA becomes the principal 
load carrier in bones of older adults. 

Reductions in estrogen and testosterones lead to in- 
creased bone turnover and an imbalance in remodeling. 
The accompanying weakening of a TA increases the sus- 
ceptibility of bone to fracture H. The management of 
osteoporosis can be greatly aided by the availability of 
characteristics that can identify this weakening. 




FIG. 1. A cross section of the TA from a 70 year old fe- 
male. Observe that the trabeculae form a disordered cubic 
network. The network loses connectivity with aging. 

Bone Mineral Density (BMD), or the effective bone 
density is the principal clinical measurement used as a 
surrogate of bone strength |1C|] . It is estimated using 
x-ray imaging via an evaluation of the opacity of bone. 
Large clinical studies have shown an exponential increase 
in fracture incidence with the loss of BMD 1 1 1 1 . Among 
other features known to be relevant for bone strength is 
the architecture of the porous bone. Structural prop- 
erties of a TA such as the average width of the tra- 



beculae, mean trabecular spacing, connectivity of the 
network and its fractal dimension have been proposed 
as additional surrogates of the ultimate (or breaking) 
strength of bone |7]|3],[12| . Alternative approaches to es- 
timate bone strength involve the use of finite clement 
computations |13| and cellular models |14|. All previous 
analyses have been carried out on TAs from bone sam- 
ples. This makes it nearly impossible to isolate effects of 
individual factors on bone decay. 

We introduce a simple mechanical model of a TA to 
complement these analyses. Since model parameters can 
be varied independently we expect the task of elucidat- 
ing the essential differences between healthy and osteo- 
porotic bone to be simplified. Possible diagnostic tools 
for osteoporosis can then be identified using their ability 
to quantify these differences. 

Images (e.g., Figure 1) suggest that disordered elas- 
tic networks p5[ may be used to model TAs. Our ini- 
tial studies are conducted on a square network of linear 
springs. An unstressed configuration is generated by dis- 
placing the vertices of the square grid (of side D) ran- 
domly by an amount less than DA. Externally imposed 
deviations increase the potential energy of the network 
via a combination of elastic (\k{Al) 2 ) and bond bending 
(^k(A9) 2 ) contributions. Here Al and A8 are changes 
in the length of a spring and the bond angle between 
adjacent springs respectively Q. The differences be- 
tween trabeculae are modeled by assigning random val- 
ues for k and k; specifically k £ [k (l — r] e ) 7 k (l + rj e )\ 
and k £ [ko(I — 77f,),Ko(l + %)], where r\ e and r\\, are 
predetermined. Since the network is constrained to lie 
on a plane, torsional effects are ignored. The following 
conditions are included to model known features of bone. 

• As discussed above, the fracture criterion for tra- 
beculae should be based on the level of strain Q • 
Following mechanical studies of fracture, we assume 
that fracture strains 7 are distributed on a (two- 
parameter) Weibull distribution llqjl^], with cu- 
mulative probability 



C(7) = 1-exp -(7/7 e ) 



(1) 



An elastic element which is strained beyond its 
fracture strain is removed from the network along 
with all bond-bending forces it contributes to. The 
bonds are assumed to fracture when changed be- 
yond an amount distributed on a second Weibull 
distribution (with parameters 7^ and m); the bond 
and its side with the smaller spring constant are 
removed from the network. 

Osteoporosis is modeled by a random removal of 
springs from the network. The probability of re- 
moval, u, is used to quantify the "level of osteo- 
porosis;" links eliminated at a given value of v are 



not re-introduced later. The surviving springs are 
assumed to retain their strength. 

Clinical studies show a dramatic increase of bone 
strength (disproportionate to the increase of BMD) 
following therapeutic intervention |FjJ. It is ex- 
plained by assuming that bone remodeling is pref- 
erentially carried out in locations of high stress 
(Wolff's Law @). In the model, elasticity of the 
spring under maximum strain is increased by a 
small fraction p; those strained by a fraction / > /o 
of the maximum are strengthened by p(jZf) ■ 
This algorithm is implemented over repeated ap- 
plication of a given external strain. 
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FIG. 2. Stress distributions on 60 x 60 networks with 
D = 10.0, A = 0.1, k = 1, k = 5, r] b = 0.5, r) a = 0.5, 
7e = 0.05, 7b = 0.1, and m = 5. For clarity only the com- 
pressed springs are shown, and elements under higher stresses 
are represented by darker hues, (a) When a network with 
v = 15% is strained by C — 4.0, the stress backbone is 
evenly distributed, (b) In contrast, the stress backbone when 
v = 30% (( = 4.0) consists of a few coherent pathways. The 
crosses denote locations with large fractures in the x-direction. 

Below, we discuss properties of networks subjected to 
uniform strain £, reached through a sequence of small, 
adiabatic increments |1S|] , When elastic elements are re- 
moved from the network due to fracture, equilibrium is 
recalculated prior to increasing £. The sides of the net- 
work are constrained to move vertically, as is the case for 
a TA because it is connected to the cortex. The force 
T(£) required to sustain a given strain is estimated by 
adding the vertical forces on the upper surface. Results 



analogous to those presented are observed under four- 
point bending which is a leading mechanism of fracture 
of long bones J2(| ■ 

The stress-strain relationship T(£), which is initially 
linear becomes nonlinear beyond the "yield point" and 
reaches a maximum (ultimate stress) T max prior to fail- 
ure of the network. The yield point coincides with the 
first fracture of springs. Indirect evidence has been pre- 
sented to suggest that yielding of a bone coincides with 
cracking of trabeculae [El| . Failure of a network is accom- 
panied by a crack propagating across the entire network. 

Figure g(a) shows the stress distributions on a 
"healthy" network strained below yield. Springs expe- 
riencing large stresses (forming the "stress backbone") 
are distributed evenly throughout the network. The 
histogram of stresses for this configuration contains 
a broad peak, similar to histograms for elastic net- 
works J22). When Q is increased, the stress backbone 
occupies a smaller subset following fracture of a group 
of springs. Similar changes have been observed in elas- 
tic networks |23fl and in finite element computations on 
digitized images of bone |U| . 

The nature of the stress backbone is very different 
in weaker networks; i.e., those with larger values of v. 
As shown in Figure §(b), the stress backbone for such 
networks consists of a few coherent pathways. Beyond 
v > Vq rs 20% there is no peak in the stress-histogram. 
We expect these differences to prove useful in identifying 
new diagnostic tools for osteoporosis. 

In Figure g, the X's denote long horizontal fractures; 
specifically, locations where four or more consecutive ver- 
tical bonds are absent. They prevent the propagation 
of stress in significant vertical slices, even when bond- 
bending forces are included. When v is large, the stress 
backbone can be seen to avoid these regions, and conse- 
quently contains only a small fraction of available bonds. 

The last observation can be used to estimate the decay 
of bone strength with increasing v. Since an externally 
applied stress passes through every horizontal layer, we 
need to consider stress propagation on each one dimen- 
sional slice. This "chain-of-bundles" model O] further 
assumes that the stress on each edge of a fracture con- 
taining k sites is enhanced by a factor (1 + -|); i.e., the 
edges share the load assigned to bonds in the fracture. 

Consider an elastic network of size N x N. The prob- 
ability of forming a fracture of length k (with links on 
either side) is v k {\ — v) 2 . Hence the size of the largest 
fracture k m (which occurs with a probability of ~ 1 24 1) 
can be estimated by 



N 2 v k ™{l-vY 



1. 



(2) 



For sufficiently large fractures k m ~ —2 In N/ In v. In the 
chain-of-bundles model, the stress on each side of this 
fracture is (1 + %p- )■£, where T is the externally applied 
stress. When this stress is the typical breaking stress 



of an elastic element, the fracture will propagate. Thus 
Tmax is related to v through 



T_~ '-- In* 



(3) 



Figure |3| shows that the decay of T max in model networks 
is consistent with this expression. 
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FIG. 3. The mean and standard error of the ultimate stress 
T max for five networks with identical control parameters as a 
function of v. The dashed line is the least squares fit. 

A similar analysis of "penny-shaped" fractures in three 
dimensional networks gives 



T„ 



V^lrTI 



(4) 



Results from mechanical studies of bone samples (with 
an assumed power law relationship) show that T max ~ 
(BMD)~ a , where a w 2.6 J25|. There is no theoretical 
basis for this power law, and Eqn. (0) also provides a 
satisfactory fit to the data. More sensitive experiments 
need to be conducted to discriminate between these pos- 
sibilities. 

We finally discuss "therapeutic regeneration" (intro- 
duced above) of a network. There has not been quanti- 
tative studies on this issue, but it has been suggested that 
therapeutic regeneration will be more effective in increas- 
ing the strength of healthy bone f25|j . This suggestion is 
based on an assumed absence of hysteresis with changes 
in BMD. The validity of this proposal can be tested in 
the model system [g6| . Figure |] shows the fractional in- 
creases of T max experienced by two networks that are 
formed by degrading the same lattice by v = 10% and 
v = 20%. As seen in clinical studies, there is a dramatic 
enhancement of bone strength (~ factor 2), dispropor- 
tionate to the increase of BMD (~ 2%). For the exam- 
ple shown, the weaker network is seen to strengthen by a 
larger factor under increase of BMD. There are large fluc- 
tuations of the enhancement between distinct networks, 
and in occasional examples the growth of the stronger 
network is larger. 

In this Letter, we have argued that an extension of 
elastic networks can be used to model the inner porous 



regions of bone. Since the latter is the principal load car- 
rier in bones of older adults, the model can prove useful 
in identifying characteristics of osteoporotic bone. Eqn. 
(^) is compatible with published data, but more sensitive 
tests are needed to discriminate it from a possible power 
law decay. Examples from the model provide evidence 
that bone regeneration under Wolff's law is not compat- 
ible with the absence of hysteresis in bone strength with 
bone density. In fact, generally, regeneration is more ef- 
ficient in strengthening networks that are weaker. 
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FIG. 4. The fractional enhancement of T ma x as a function 
of the percentage increase of BMD under therapeutic regen- 
eration (with /3 — 2, fo — 0.8) for two networks with v — 10% 
(dashed line) and v = 20% (solid line). 

We propose to use the model identify new diagnos- 
tic tools for osteoporosis, using characteristics that can 
differentiate between stress backbones (Figure ||). Pre- 
liminary studies indicate that the ratio of responses of a 
network to static and periodic strain is a suitable surro- 
gate of bone strength |27]] . 
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